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Abstract 

We construct models of many-particle quantum graphs with singular two-particle 
contact interactions, which can be either hardcore- or 5-inter actions. Self-adjoint 
realisations of the two-particle Laplacian including such interactions are obtained via 
their associated quadratic forms. We prove discreteness of spectra as well as Weyl 
laws for the asymptotic eigenvalue counts. These constructions are first performed 
for two distinguishable particles and then for two identicle bosons. Furthermore, we 
extend the models to N bosons with two-particle interactions, thus implementing the 
Lieb-Liniger model on a graph. 
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1 Introduction 



This paper is the second in a serious of papers developing models of interacting, non- 
relativisitc many-particle systems on (compact) graphs. The underlying one-particle quan- 
tum graphs are well established models for a variety of problems in quantum mechanics. 
Interest in such one-particle models surged after Kottos and Smilansky [KS99bJ had dis- 
covered that the eigenvalues of quantum Hamiltonians describing single particles on graphs 
possess the same correlations as eigenvalues of random hermitian matrices. Henceforth, 
quantum graphs proved to be very successful models in the area of quantum chaos [GS06J 
and beyond |EKK+08j . 

In our first paper [BK11J we introduced systems of two particles on compact metric 
graphs. We identified self-adjoint realisations of the two-particle Laplacian that describe 
singular two-particle interactions that are located in the vertices of the graph. In that 
context we analysed the Laplacians indirectly, by first constructing closed, semi-bounded 
quadratic forms and then identifying the self-adjoint operators that are uniquely associated 
with the quadratic forms. Using quadratic forms allowed us, furthermore, to prove that the 
Laplacians have compact resolvents and thus possess purely discrete spectra. Moreover, 
with a bracketing argument we were able to prove a Weyl law for the asymptotic eigenvalue 
count. 

The goal of this paper now is to introduce two-particle interactions of a different kind, 
modelling short-range interactions in terms of singular contact interactions that are for- 
mally given by a Hamiltonian 



Here x and y are coordinates for the positions of the two particles on the edges of the 
graph and a 6 R is a coupling parameter, such that a > corresponds to repulsive 
interactions. Contact interactions of this kind play a prominent role in the description of 
Bose-Einstein condensates (BEC), leading to Gross-Pitaevskii equations, in plasma physics 
and in solid state physics (see, e.g., |LSSY05[ lCCG + li] ). They are also of interest from 
a mathematical point of view: models of (an arbitrary number of) particles on the real 
line with Hamiltonian ( 11. ip are completely solvable in the sense that the eigenfunctions 
can be constructed explicitly, see |Yan67l IAFK02t IAFK04] . As in the one-particle case the 
connectivity of a graph, however, adds sufficient complexity to the problem to render it not 
solvable in that sense. Therefore, bosonic many-particle systems on graphs are expected 
to show generic many-particle properties of systems confined to one spatial dimension. 

In this paper we shall follow our previous approach [BK11] in that we first introduce 
suitable quadratic forms, show that they are closed and semi-bounded, and then iden- 
tify the corresponding self- adjoint operators. These operators shall always be two-particle 
Lapacians, however, with suitable domains that implement singular two-particle contact 
interactions. The resulting operators are then rigorous versions of the formal Hamilto- 
nian (II. ip . In order to achieve this the domains of the operators are required to contain 
jump conditions along diagonals x = y for derivatives of the functions in the domain. A 
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similar construction was already given by Harmer |Har07t IHar08j , who considered two par- 
ticles on a star graph (with semi-infinite edges). Our construction works for any compact 
metric graph and allows for straight-forward generalisations to somewhat broader classes 
of singular contact interactions of which the 5-type interactions (11 .ip form a prominent 
subclass. 

The paper is organised as follows: In Section [2] we review important properties of one- 
particle quantum graphs and introduce relevant concepts and notations that we are using 
in the following sections. Section [3] is then devoted to the construction of the contact 
interactions via quadratic forms for two distinguishable particles on a general compact, 
metric graph. In that context we encounter the problem of elliptic regularity |Dob05] in 
the same way as previously |BKllj . suggesting an analogous solution for a similar class of 
boundary conditions. In Section H] we first implement a bosonic realisation of a particle 
exchange symmetry and hence obtain a rigorous construction of the bosonic version of the 
operator (II. ip . We then extend this construction to N bosons on a general compact graph, 
for which the equivalent to the formal operator ( II. ip is 

N 

H = - A N + a ^2 K x i ~ x i) ■ (!-2) 

i<j 

Here x\, . . . ,xn are the particle positions and — An is the Laplacian in these N variables. 
Hence, the model we are investigating is an extension to quantum graphs of the well-known 
Lieb-Liniger model |LL63j that has been studied extensively in the context of BEC. 

2 Preliminaries 

Before introducing many-particle systems on graphs we briefly describe one-particle quan- 
tum graphs. They form the basis for the tensor-product construction of many-particle 
quantum systems on graphs. 

The classical configuration space of a quantum graph is a compact metric graph, i.e., 
a finite graph T = (V, £) with vertices V = {v i, . . . , vy} and edges £ = {ei, . . . , eg}. The 
latter are identified with intervals [0,Z e ], e = 1,...,E, thus introducing a metric on the 
graph, see |KS99bl IKS99al IKuc04l CSP6] for details. 

Functions F on the graph are collections of functions on the edges, 

F — (A) ■■■■> Ie) , with f e :[0,l e )^C, (2.1) 

so that spaces of functions on T are (finite) direct sums of the respective spaces of functions 
on the edges. The most relevant space is the one-particle Hilbert space 

E 

^ 1 = L 2 (r):=0L 2 (O,/ e ), (2.2) 

e=l 

and all other spaces are constructed in a similar way. 
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The one-particle Hamiltonian is a Laplacian, acting as a second-order differentiation, 

- AiF = (-f[ f , ...,-/*), (2.3) 

on F G C°°(r). We here use the index to indicate that this is a one-particle Laplacian. 

Domains of self-adjoint realisations of the Laplacian are characterised in terms of bound- 
ary conditions in the vertices. These require boundary values 

F bv := (/i(0), . . . , MO), Mh), . . . , f E (l E )) T G C 2E , (2.4) 

of functions and (inward) derivatives, 

FL := (f[(0), f' E (0), -fM, -f' E ilE)) T 6 C 2E , (2.5) 

as well as projectors P and Q = \ie ~ P acting on the space C 2E of boundary values, and 
self-adjoint endomorphisms L of ranQ C C 2E . 

The self-adjoint realisations of — Ai can be identified via the quadratic forms that are 
uniquely associated with them |Kuc04] . 

Theorem 2.1 (Kuchment). The quadratic form 

Qp!l[F] = [ |V/| dx - (F bV} LF bv } C 2 E = V f E \f e (x)\ 2 dx - (F bvi LF^)^ , (2.6) 
Jr Jo 



e=l 



with domain 

V QW = {F eH\T)- PF bv = 0} (2.7) 

is closed and bounded from below. The unique, self-adjoint and semi-bounded operator 
associated with this form is the one-particle Laplacian — A\ with domain 

Pi(P, L) = {F G H 2 (T); PF bv = and QF' hv + LQF bv = 0} . (2.8) 

A two-particle quantum system requires the tensor product of two one-particle Hilbert 

spaces, 

U2:=Ui®Ui . (2.9) 
For a quantum graph this means that 

E E 

H 2 = (0 L 2 (0, i e )) ® (0 L 2 (0, U) , (2.10) 

e=l e=l 

such that vectors ^ G T-L% are collections \1/ = (ip ei e 2 ) of -E 2 functions defined on the 
rectangles -D eie2 = (0,/ ei ) x (0,/ e2 ). Their disjoint union is denoted as 

D T = \jD eie2 , (2.11) 

eie 2 



so that one may view %2 as 



L 2 (D r ):=Q)L 2 (D eie2 ). 



(2.12) 



We shall use a similar notation for other function spaces. 

As a differential operator, the two-particle Laplacian acts as 



(-A a tt) 



d 2 



' l Pe 1 e2 

dx\ 



& 2 



(2.13) 



dx 2 



and hence has the same form as a Laplacian in M 2 . Defined on the domain C^°(Dr), this 
operator is symmetric, but not self-adjoint. 

Self-adjoint realisations of the two-particle Laplacians can either represent non-inter- 
acting particles, or introduce two-particle interactions via boundary conditions. A par- 
ticular class of singular two-particle interactions that are localised in the vertices was 
established in |BK11] . Here we shall introduce two-particle contact interactions that are 
localised on the edges. 

3 Contact interactions on a general compact graph 

The interactions we have in mind are intended to model two point-like particles on a graph 
that interact when they hit each other, i.e., when they are located in the same position. 
This requires, in particular, that they are on the same edge. Hence, the subset of the two- 
particle configuration space Dr (12. lip where these interactions take place consists of the 
diagonals of the squares D ee . Singular interactions require a dissection of the configuration 
space along these subspaces, and suitable matching conditions for functions and their 
derivatives along the boundaries introduced by the dissection. We therefore define the 
'dissected' configuration space 



where Z)+ = {(x, y) G D ee ; x > y} and D ee = {(x, y) G D ee ; x < y}. Functions on are 
denoted as \I/ = (i/; eie2 ). The components ip ei e 2 for t\ ^ e<i are defined on D eie2 , whereas 

i^ee = (ipte^ee) with V>« defined on Df 

The two-particle Hilbert space "H 2 (12. 10)) can then also be viewed as 






(3.2) 



Boundary values of functions \I/ G if 1 (Dp) are encoded in vectors 

*&«(3/) = {^e ie2 ,bv{y)) and %v(y) = (^^(y)) 



(3.3) 
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We distinguish components with e\ 7^ e 2 from those with e\ = e2, as in the latter case ad- 
ditional boundary values along diagonals have to be taken into account. More specifically, 
when ei 7^ e 2 we define 



ip ei e 2 ,bv(y) ■-- 



l ei ip eie2 (l ei y,0) 



/ 



and 



^' ei e 2 ,bv(y) ■-- 



^ ei e 2 ,x(0,/ e2 y) \ 

/l e2 ijj eie2tX (i ei , le 2 y) 

_ lei'll , e 1 e2,y(le i y?Q) 
162,?/ 

(l ei y,l e2 )J 

(3.4) 

where y G [0, 1]. When e\ = e 2 boundary values along the diagonals of the squares D ee have 
to be added, including those for derivatives. Noting that the inward normal derivatives 
along the 'diagonal' part of the boundary of Df e are 



± 

ee,n 



^2 (Vee,a; ^ee,y) 



(3.5) 



we set 



1pee,bv(y) 



( Vfe^-(0,l e y) \ 

y/kip+(l e ,l e y) 

^te(hy,o) 

h^eeihyje) 

htpte( l ey,hy) 

\Vh^ee( l ey, Uy) J 



and 



/ Vkll>ee,x( Q , l ey) \ 
-VhlPte,x( l eil e y) 

Vh^te,y( l ey,°) 
-Vh^ee,y(hyJe) 

y^h^nihyjey) 

\s/W e 'il)- {l e y, l e y)J 



(3.6) 



for y G [0, 1]. Altogether, the vectors ( 13. 3 p of boundary values have n(E) := 4E 2 + 2E 
components. 

As a next step we introduce the bounded and measurable maps P, L : [0, 1] — > M(n(E), C) 
that are required to fulfil 

1. P(y) is an orthogonal projector, 

2. L(y) is a self-adjoint endomorphism on kerP(y), 

for a.e. y G [0, 1]. We then introduce two bounded and self-adjoint operators, II and A, on 
L 2 (0, 1) ®C n{E) . They are defined to act as (ILx)(v) '■= p (y)x(y) and (Ax)(y) ■= L{y)x{y) 
on x G L 2 (0,1) ®C< E \ 

Our aim is to obtain self-adjoint realisations of the two-particle Laplacian — A 2 , see 
f)2.13p . that represent two-particle contact interactions and are extensions of — A 2>0 defined 
on the domain C^(D^). In analogy to the one-particle case (12.81) . as well as the case of 
singular interactions covered in [ BKllj . their domains should be given in the form 



V 2 (P,L) := G H\D*); P(y)V bv (y) = and 

Q(y)% v (y) + L{y)Q{y)m hv {y) = for a.e. y G [0, 1}} 



(3.7) 



In order to clearly distinguish the boundary conditions that induce contact interactions 
from other kinds of boundary conditions we rearrange the order of terms in the boundary 
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vectors. We first list, for each edge e, the lower two boundary values in ( 13. 6p . and then, 
for each pair (ei,e 2 ), either the four components in (13 .4p or the upper four components of 
(13. 6p . That way one achieves a decomposition of the space of boundary values according 
to 

C ^ ^ = Vcontact © Vvertex ■ (3-8) 

Here V con tacti with dimV cont act = 2E, contains the boundary values (I3.6P along diagonals, 
whereas V vertex , with dim V vertex = 4E 2 , contains the remaining boundary values (13.4)) and 
(13. 6p . which are associated with vertices. A separation of contact interactions from any 
other boundary effects requires to choose P and L as block-diagonal with respect to the 
decomposition (I3.8p . From now on we assume this to be the case. 

For the restriction of P and L to V ver t ex we assume the same conditions as in |BKllj . For 
most purposes, however, it is sufficient to suppose that there are no two-particle interactions 
in the vertices. In [BKllj the non-interacting boundary conditions were characterised as 
follows: The restrictions of P and L to V vertex are independent of y and block-diagonal 
with respect to a decomposition of V verte x according to the index e 2 in (13.31) . 

The restriction of P and L to V conta ct should, first of all, be block-diagonal with respect 
to a decomposition of that space according to the edges in order to avoid 'contact' interac- 
tions across edges. Further restriction are not necessary, but we identify the following two 
cases as of particular interest because they correspond to a Hamiltonian of the form (11.11) . 

Definition 3.1. Let a : [0, 1] — > M be Lipschitz continuous. Then a contact interaction is 
said to be of 

(i) 5-type with (variable) strength a, if \& G H 2 (Dy) is continuous across diagonals, 

TpteihV, hy) = 1p7e( l eV, hy) , (3-9) 

and satisfies jump conditions for the normal derivatives, 

*Pte,n(ley, hy) + C e ,n(^2/, hy) = -^=a(y)i)f e (l e y , hy) , (3.10) 

(ii) hardcore type, if it satisfies Dirichlet boundary conditions along diagonals. 

We remark that contact interactions of the 5-type can be seen as a rigorous realisation 
of a Hamiltonian 

- A 2 + a(y) 5(x - y) . (3.11) 

The case a(y) > for all y G [0, 1] corresponds to repulsive interactions and is the most 
relevant case for models of actual particles on a graph. Hardcore interactions follow from 
such a formal Hamiltonian in the limit a — > oo. 

Following our intention to represent domains of two-particle Laplacians in the form 
( 13. 7p we have to choose the maps P and L in such a way that their restrictions to the 
edge-e subspace of V contac t are 

1/1 -1 



con I ad j \ ../ / ^ 



W = 2 -1 f (312) 
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and 

L C ontact,e{y) = ~^ 01 ^ 1<2 ' ( 3 ' 13 ) 

in order to generate 5-type contact interactions. Hardcore interactions require the choice 

P contact — 1 and L contact = 0. 

Our approach to self-adjoint realisations of the Laplacian uses suitable quadratic forms, 
which are uniquely associated with these operators. 

Proposition 3.2. Assume that the maps P,L : [0,1] — > M(n(E),C) are bounded and 
measurable. Then the quadratic form 

Q i p ) M = (W,W) L 2 {D * r) - [ (*to(y),L(y)*Uv))c«*)dy , (3.14) 

Jo 

with domain 

V Q m = {V G H\D*); P^^y) = for a.e. y G [0, 1]} (3.15) 
is closed and semi-bounded. 

Proof. The proof follows by using the same steps as in the corresponding proof in |BK11] . 
The only consideration that has to be added concerns the upper bound 

i 

(Vbv(y),L(y)Vbo(y)) C n(E) dy < L max \\^ bv \\ 2 L2{ ^ mcn(E) . (3.16) 
To estimate the right-hand side one requires the bound, 

II^IIW)^) < K Wv{d$ + 5 UV*||£ 2(Df) ) , (3.17) 

to hold for all 5 < 5q, where K, 5q > 0. The contribution from the rectangles D ei£2 (with 
e! / e 2 ) in the decomposition ( 13. ip can be dealt with as in |BK11] and is based on a result 
in |Kuc04j . For the triangles Df e we note that close to the corners with angles tt/A this 
method fails. However, one can always reflect functions i/jf e across edges, define them on 
suitable squares and then apply the bound as before for the rectangles. The proof then 
continues as in [BK11] . □ 

According to the representation theorem for quadratic forms (see, e.g., |Kat66j ) there 
exists a unique self-adjoint and semibounded operator H with domain V(H) C Vqw that 

is associated with the quadratic form Q PL . It is not immediately clear, however, that the 
functions in T>(H) possess if 2 -regularity. If this is the case we say, for short, that the 
quadratic form is regular. We note that a self-adjoint realisations of — A2 with domain 
(13. 7p would correspond to a regular form. 

Under an additional (mild) assumption a regular quadratic form indeed leads to a 
two-particle Laplacian with domain (13.71) . 
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Proposition 3.3. Suppose that the map P is of class C 1 and that the quadratic form Q PL 
with domain Vn(2) is regular. Then the unique, self-adjoint and semibounded operator that 
is associated with this form is the two-particle Laplacian — A 2 with domain V 2 (P,L). 

Proof. The proof can essentially be taken over verbatim from the corresponding proof in 
|BKllj . It is based on the representation theorem for quadratic forms, which implies that 
for each \1/ G T)(H) there exists a unique \ G L 2 (D^) such that 

Qg ) L [vl> ! $] = ( x ,$), V0GP Q(2) . (3.18) 

When $ G C^(D^), an integration by parts of (I3.18P implies that H acts as a two-particle 
Laplacian — A 2 . In the general case of a \& G T)q(2) the integration by parts yields an 
additional boundary term, 

-/ (*'bv(y) + L{y)Vbv(y),®Uy))c"W d y , (3.19) 

Jo 

that is required to vanish. Following Lemma 3.13 in |BKllj . which has an immediate 
generalisation to the present case, the set $ G T>q<2)} is dense in ker II C L 2 (0, i) n ( E '. 

Hence, ^ + A^ bv G ker IL 1 , or 

Q(y)%M + Q(y)L(y)^ bv (y) = ■ (3.20) 
This condition finally implies that T>(H) = T>2(P, L). □ 

As mentioned above, the quadratic forms in Proposition 13.21 are not necessarily regular. 
Since our focus is on contact interactions of 5- or hardcore- type, it is sufficient to consider 
these cases. However, as in |BK11] we have to add one additional assumption on the 
restrictions P ver t of the projectors P to V verte x- Splitting V veriex into the two subspaces 
spanned by the upper two and the lower two components of (13. 4p as well as the upper two 
and middle two components of f )3.6p . respectively, we require P ver t to be block-diagonal 
with respect to this decomposition. 

This then leads to the main result of this section. 

Theorem 3.4. In addition to the assumption made for the maps P and L above, suppose 
that P is of class C 3 and L is Lipschitz continuous. Furthermore, for y G [0, ei] U [I — 62, 1] 
with some ei,e 2 > assume that the restriction of P to V ver t ex is diagonal with diagonal 
entries zero or one and, in the case of 5 -type interactions, that a(y) = «o > for those y. 
Then the quadratic form Q^ L is regular. 

Proof. First note that it is enough to show regularity near the corners of D ee = D+ e U D~ e 
adjacent to the diagonal. The regularity away from the diagonal of D ee as well as regularity 
in the rectangles D eie2 with e\ 7^ e 2 was already established in |BKllj . In addition, the 
regularity along the diagonals in the interior of D ee can be readily established using the 
same methods as in |BKllj . 
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The assumptions made on P imply that on the edges of the squares D ee the functions 
in T>q(2) satisfy either Dirichlet- or Neumann boundary conditions near the corners. Along 
diagonals we consider the projections 

I (3-21) 

tpee,F(x, y) : = - [lpee{x, y) ~ 1pee{y, x)] . 

The goal is to show that, close to the corners, both i/) ee> B and ip ee ,F are of class H 2 . For that 
purpose one introduces suitable cut-offs that restrict the functions (I3.2ip to neighbourhoods 
of the corners. This eventually implies that tp ee G H 2 (D* ee ). 
We recall the conditions (13.101) which imply that 

dn^ e , B ~ ^=<b = , (3.22) 

on the diagonal. Hence, ipf e B satisfies (variable) Robin boundary conditions on the di- 
agonal. By construction, i/jf e F vanishes on the diagonal so that near the corners of Df e 
adjacent to the diagonal, where a is supposed to be constant, ^ eeB / F satisfies a combi- 
nation of Dirichlet-, Neumann- or standard Robin-boundary conditions. In all such cases 
regularity is well known to hold |Dau88[ IGri85j . □ 

One naturally expects the two-particle operators representing contact interactions to 
possess purely discrete spectra of the form Ao < Ai < A2 < . . . (i.e., eigenvalues are 
counted with their multiplicities and do not accumulate at any finite value). Moreover, 
the asymptotic distribution of eigenvalues, as given by the asymptotic behaviour of the 
eigenvalue-counting function 

N{\) := A n < A} , (3.23) 

should follow a Weyl law. We shall now prove a Weyl law for repulsive contact interactions. 
This includes hardcore- and ^-interactions with a > 0. The general requirement is that 
L contact is negative definite (compare (I3.13P ). 

Proposition 3.5. Let (—A2,T>2(P,L)) be a self-adjoint realisation of the two-particle 
Laplacian with repulsive contact interaction as described in Proposition \3.3i Then this 
operator has compact resolvent. In particular, its spectrum is purely discrete and only 
accumulates at infinity. Furthermore, the counting function (I3.23P obeys the Weyl law 

C 2 

iV(A) ~ — A, A ^00, (3.24) 

47T 

where £ = J2e=i b is the total length of the graph. 
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Proof. The Hilbert space H 1 (D^) is compactly embedded in L 2 (D^). Since the form norm 
|| • ||g(2) is equivalent to the H 1 (D^)-norm, the Hilbert space (XV>(2), || • ||q(2)) is also com- 
pactly embedded in L 2 (D^). Hence the operator associated with the quadratic form has 
compact resolvent [Dob05j . 

The Weyl law follows from a standard bracketing argument [RS78J based on a compar- 
ison with two suitable operators (quadratic forms), see also [BE09, BKllJ. 

The first operator, (— A 2 , V 2 (P D , Ld)), is the Dirichlet-Laplacian, and is characterised 
by the projector P D = 1 as well as L D = 0. Given an operator (— A 2 , T> 2 (P, L)) the second 
comparison operator, (— A 2 , T>2(Pr, Lr)), is the Laplacian given by the projector Pr = 
as well as 

L R = diag(A 1 _^A,0 1 _^0) , (3.25) 

4E 2 — times 2E— times 

where A = ||A|| op . The associated quadratic forms therefore satisfy the following inclusions 
of their domains, 

V 2 {P D , L D ) C V 2 (P, L) C V 2 (P R , L R ) . (3.26) 
Hence |RS78j . it follows that the related eigenvalue-counting functions satisfy 

N D (X) < N(X) < N R (X) . (3.27) 

As both No and Nr satisfy the Weyl law (13.241) . the same asymptotics hold for N(A). □ 



4 Contact interactions for bosons 

So far we assumed to have two non-identical particles on a graph. We now implement an 
exchange symmetry for two identical bosonic particles. Their states are described in the 
symmetric two-particle Hilbert-space H 2t b — Hi ®s H\. The orthogonal projection n# 
from 7-L 2 = L 2 (D^) to the bosonic subspace 

L 2 B (D*) = U B L 2 (D*) (4.1) 

acts on components of \I/ = (VWea) with e\ ^ e 2 as 

eieaC^en Ve-i) — „ yPe\ei{^e.\i Ve-i) tyeiex {lle2J -^ei)) j (4-2) 

whereas on the components with e± = e 2 its action reads 

(U B ^)f e (x e ,y e ) = ^±(x e ,y e ) + ^ e (y e ,x e )) . (4.3) 

Due to this symmetry it would be sufficient to keep only components with e\ < e 2 in 
addition to the diagonal components with e\ = e 2 . For simplicity, when comparing to the 
previous section we, however, keep all components. We then denote the images of function 
spaces under the projection to their bosonic subspaces as, e.g., 

H%{D*) = H m {D* T ) n H 2 , B . (4.4) 
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We note that whenever \1/ G H2,b is i n H 2 (Dy), the underlying symmetry implies the 
relations 

^eie2,x(^ei > ?/e 2 ) = "0e2ei,j/(2/e2) -^ei) and 1peie2,xx(%e 1 , ye 2 ) = V^ei ,yy (Z/e 2 ? ^ej ) j (4-5) 

when ei 7^ e2, as well as 

V^zOe, y B ) = ^yiVe, Ze) and ll>%,JXe, y e ) = 1pZ,yy(Ve, *e) • (4-6) 

Due to the bosonic symmetry it is possible to reduce the number of components in the vec- 
tors of boundary values (13. 3p . When ej 7^ it suffices to keep the upper two components 
in each of the vectors (13.41) . whereas for e\ = t2 we use 



^ee,bv(y)= I VklPteiteihy) and ^^(y) = -^/%1p+ x (l ej l e y) | . (4.7) 

^tei^yjey)/ \y/W e ll)+ (l e y,l e y) 



with y G [0, 1]. The space of boundary values therefore has dimension ub(E) = 2E 2 + E, 
and decomposes in analogy to (I3.8p . 

We also need the bounded and measurable maps P, L : [0, 1] — >■ M(n B (E), C), where 

1. P(y) is an orthogonal projector, 

2. L(y) is a self-adjoint endomorphism on kerP(y), 

for a.e. y G [0, 1]. The space C ns ^ of boundary values decomposes in the same way as 
(13. 8p . however, the edge-e subspaces are now one-dimensional. This forces the equivalent 
of (I3.12p to be a projector on C and to take values 



contact, e 



(y)e{0,l}. (4.{ 



Here P con tact,e(y) = 1 corresponds to a Dirichlet condition in the point (l e y,l e y) along the 
diagonal, whereas P CO ntact,e(y) = imposes no condition. Hence, when 5-type interac- 
tions are considered we choose P CO ntact,e(y) = 0, and in the case of hardcore-interactions 
P C ontact,e(y) = 1 is chosen. Likewise, the equivalent of (13.131) is 

L C ontact,c{y) = ~~ Oc{y) (4.9) 

for (^-interactions, and L conta ct,e(y) = for interactions of hardcore type. 
We can now set up the following quadratic form, 

Q i p\ B VP]=2(m x ^ x ) L 2 {D , ) -2 [ (* to (j/), J L(y)* to (y))c» fl( i S ) dy , (4.10) 

Jo 

with domain 

V Qm , B = {V G Hl(D*); P(y)V bv (y) = for a.e. y G [0, 1]} . (4.11) 
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As this is the restriction of a quadratic form on L 2 (.Df) to L 2 B (D^), all results of Section [3] 
carry over: Propositions 13.21 and 13.31 imply that the quadratic form is closed and semi- 
bounded; when P is of class C 1 and the form is regular, the associated self-adjoint operator 
is the bosonic two-particle Laplacian — A2,_b with domain 

T\b(P, L) := {* G H 2 B (D* r ); P(y)^ bv (y) = and 

QiyWM + L(y)Q(y)V bv (y) = for a.e. y G [0, 1]} . 1 ' ' 

According to Theorem I3.4[ when the consitions of that theorem are fulfiled the quadratic 
forms leading to 5-type and hardcore-interactions are regular. 

We remark that for ^-interactions one can use the decomposition ( 13. 8 j) of the space of 
boundary values and the explicit expression (14. 9p to rewrite the quadratic form as 

Qp ] L B bP} = 2 (^x, ^x)l%(d*) ~ 2 / (^bv,vert(y), L vert {y)^ bv>vert {y)) c2E 2 dy 

E x (4.13) 

+ ^2 a ^ \V^e1pte( l ey,l e y)\ 2 dy . 
e=l J ° 

In the same way, the form domain takes the form 

V qW ,b = G H B (D*); P ver t{y)^bv,vert{y) = for a.e. y G [0, 1]} . (4.14) 

Due to the bosonic projection n#, which commutes with any of the two-particle Lapla- 
cians, asymptotically half of the spectrum of a Laplacian is projected to the bosonic Hilbert 
space Ti<2,Bi so that Proposition 13.51 implies the Weyl law 

jV B (A)~^A, A^oc, (4.15) 

87T 

for the asymptotics of the eigenvalue count restricted to %2,b- 

Our goal now is to study bosonic many-particle systems on graphs. Eventually, these 
have to be described in the bosonic Fock space over the one-particle Hilbert space. Since 
it suffices, however, to consider each iV-particle space separately, we here only consider a 
fixed particle number N. In that context we shall introduce two-particle interactions that 
are (formally) of the type, 

H N = -A N + J2 a ( x i) ■ ( 4 - 16 ) 

i<j 

Due to the bosonic symmetry, on suitable functions the quadratic form associated with 
such an operator will be 



[%H N #) HatB = -A N #) HatB 

N(N — 1) ^ fa fa 2 
+ / / ••• / o:{x e2 )\yj e2e2s _ eN [x e2 , x e2 , . . . , x eN )\ dx eN ...dx e2 . 

„. Jo Jo 

(4.17) 



e 2 ...e N 
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From (I4.16P and (14.1 7\i one concludes that contact interactions involve boundary values 
along hyper surf aces that are characterised by the fact that two particles are at the same 
position. 

The configuration space of iV (distinguishable) particles is 

D? = |J D eie2 ... eN , (4.18) 

eie 2 ...e ]v 

where -D eie2 ... ejv = (0, l ei ) x • • • x (0, l eN ). We stress that this notation includes cases where 
several particles are on the same edge, in which case the same edge appears repeatedly. The 
hyperplanes that determine contact interactions are characterised by equations x l e = x J e , 
meaning that particles i and j sit on the same position on edge e. In analogy to (13. ip . in 
order to implement contact interactions we have to decompose further; this involves 
all hyperrectangles D eie2 .„ eN that are composed of at least two coinciding edges. 

Now assume that (ni, . . . , He) is a partition of N such that there are n e particles on edge 
e. Let cr G Sn assigns labels to the N particles in such a way that cr(l), . . . , a{n e ) label the 
particles on edge e, with coordinates Xe^\ ■ ■ ■ , ie'" e ' '. Permutations tt G S ni x • • • x S nE C 
Sn of particle labels then leave the assignment to edges untouched, and there exists such 
a permutation with 

a£Wi)) < . . . < a^ (<T( " e)) , Ve G {1, . . . , E} . (4.19) 

These relations define a polyhedral subdomain of D eu „ eN . Every other permutation tt G 
S ni x • • ■ x S nE will produce a copy of that polyhedral subdomain that emerges through 
reflections in a succession of boundary hyperplanes. We will enumerate these ni!... rig! 
subdomains as with 1 < r\ < ni!. ..rig!. 

In analogy to (13.11) we now introduce the dissected hyperrectangles as the disjoint union 

K...e N = (J D l-e N ■ (4- 20 ) 
V 

The iV-particle Hilbert space Hn for (distinguishable) particles with contact interactions 
can then be defined as 

L\D^)= L 2 (£> e * ie2 ... e J . (4.21) 

eie2...ejv 

The corresponding Sobolev spaces are defined in the same way. Note that on the right- 
hand side D* ie2 ejv = Z) eie2 ... ejv when all edges in the definition of Z) eie2 ... ejv are distinct, 
i.e., when no two particles are on the same edge. With this proviso, we denote functions 
on D^ ie2 by ip ei ...e N1 which themselves consist of ni!. .. rig! components defined on the 
subdomains compare (13. ip and below. 

The projection from (14.211) to the bosonic iV-particle Hilbert space L 2 B (D^*) then 
is given by 

(n B *)e,.. e * = ^ E ^(e 1 )... 7 r( eJ v)«S) ) -^S)) • (4-22) 

7rGSjV 
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In analogy to ( 13. 9 p and (I3.10p . two-particle interactions of a 5- type (14.161) require 
boundary values of functions \1/ G Hg(D^*) and their normal derivatives along (internal) 
boundary hyperplanes of the dissected hyperrectangles D* . In addition, boundary 
conditions at vertices have to be implemented. For those purposes the most convenient 
expression for the quadratic form is an analogue of (I4.13p . 

We first introduce the vectors of boundary values in vertices. Due to the bosonic 
symmetry these can be given in the form 

xT , („,\ _ ( V^gj • • • le N 1p ei ...e N (0, h 2 Vl, ■ ■ ■ ) le N UN-l) \ /a r>o\ 

* bv,vert\y J ~ I f, 5 , // / 7 \ I ) l^- zo J 

VV <e 2 • • • l e N iPe 1 ...e N {l'e 1 , <e 2 l/l, ■ ■ ■ j ^et^VN-l) J 

and 



, I V ^ e 2 • • • ^e N ^e J _...e N ,x\. (0; ^e 2 2/l) • • • > ^e N VN-l) . -ion 
VbwertiV) = I 77 ]—.,. (11,, I „ > I ' ( 424 ) 

where ?/ = (y 1; . . . , yjv-i) G [0, l] 7 ^ 1 . On these (vertex related) boundary values the 
bounded and measurable maps P vert ,L vert : [0, l]^ -1 — >■ M(2_E Ar ,C) shall act, which are 
required to fulfil 

1. P V ert{y) is an orthogonal projector, 

2. L vert (y) is a self-adjoint endomorphism on ker P„ eri (?/), 
for a.e. y G [0, l]^ 1 . 

Boundary values on internal hyperplanes in the dissected hyperrectangles involve com- 
ponents ^ei.-.ew with a pair of coinciding edges, e^ = e^. Due to the exchange symmetry 
we can always arrange for these edges to be ex and e 2 = e.\. Permuting a given pair (e*, ej) 
to (ex, e 2 ), however, involves a change of the associated domain to some other copy 

D v , , . This means that ib% „ Ar is replaced by lb" , , , . 

e 1 —e N rei...ejv ^ J T'e 1 e 1 ...e^ 

In analogy to (I4.13P the quadratic form we wish to set up is 

Qb[^\= N Y] / ••• / \iJe 1 ...e N ,x ei (Xe 1 ,... } X eN )\ 2 dx eN ...dx ei 

N I (^bv,verti L* ver t {y)^ bv,vert) £2E N dy (4.25) 

•/[o.ip- 1 

+ iV(iV 1_ yv I «(?/!) \y/l e2 . . lZi>e 2 e 2 ...e N (le2yiJy)\ 2 d V ■ 

1 ' „ ./fcuF- 1 



e 2 ...e N 



For convenience we here used the notation ly = (l e2 yi, h 3 V2, ■ ■ ■ , h N UN-i)- 
This form shall be defined on the domain 

V Q m = {^E H l B {D^*); P vert {y)*bv,vert{y) = for a.e. y G [0, if' 1 } . (4.26) 

Using this, we can readily establish the following statements. These are immediate general- 
isations of the corresponding statements, Propositions 13.21 and I3.3[ for two bosons (N = 2). 
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Theorem 4.1. Let the maps P vert ,L vert : [0,1]^ 1 — > M(2E , C) as well as the function 
a : [0, 1] — > C be bounded and measurable. Then: 

(i) The quadratic form defined on the domain V <jv) is closed and semi-bounded. 

(ii) If P ver t is of class C 1 and the form is regular, the associated self-adjoint operator is 
the N -particle Laplacian — An with domain 

V N>B {P,L) := {V € H 2 B (D»*); P(y)^Uy) = and 

QivWbviv) + L(y)Q(y)*Uy) = for a.e. y G [0, l]"" 1 } . 

(4.27) 

Here P = P con tact © Pvert and L = L contact © L vert refer to all boundary values. 

The proof of this Theorem is an immediate extension of the proofs of Propositions 13.21 
and 13.31 as well as of Theorem 13.41 

It is also immediately clear from the proof of Proposition 13. 51 that any of the iV-particle 
Laplacians — with repulsive contact interactions have compact resolvent and hence 
possess purely discrete spectra, accumulating only at infinity. Furthermore, the eigenvalue 
counting function (compare (I3.23P and (14.151) ) satisfies a Weyl law that follows from a 
bracketing argument in the same way as (13.241) . For the case of iV distuingishable particles 
the Weyl law is 

rN 

^>~ (4 ^ r(1 + f) ^ /2 - <" 8 > 

This follows most easily from the lower bound given by the Dirichlet Laplacian as in (I3.27p . 
The bosonic case requires to desymmetrise the spectrum with respect to particle exchange 
symmetry; hence, the bosonic counting function is reduced by a factor of j^. 
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